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£S) ' Abstract 

|- An asymmetric nature of the boson 'destruction' operator a and its 'creation' partner a* is made 

apparent by applying them to a quantum state \ip) different from the Fock state \n). We show that it is 
possible to increase (by many times or by any quantity) the mean number of quanta in the new 'photon- 
ly-^ ' subtracted' state &\ip). Moreover, for certain 'hyper-Poissonian' states the mean number of quanta in 

^s^j | the (normalized) state a\ip) can be much greater than in the 'photon-added' state a^\ip). The explanation 

of this 'paradox' is given and some examples elucidating the meaning of Mandel's g-parameter and the 
£NJ ' exponential phase operators are considered. 
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CO ■ 1 Introduction 

O ; 

The non-Hermitian bosonic operators a and a^ of the harmonic oscillator, satisfying the canonical commuta- 
| tion relation [a, a* 1 ] = 1, are usually called 'annihilation' [jlj (or 'destruction' j^, ||) and 'creation' operators 

f^*) ■ (perhaps, only P A M Dirac, in his book (J], did not use these terms, introducing instead of a and a 1 the 

'complex dynamical variables' fj and ?y.) This is due to their action on the Fock (number) state \n) (eigenstate 

of the number operator n — a) a) Q 
±\ ' a\n) = V^\n-l), (1) 

a f |n) =Vn+l \n + l), (2) 

Therefore, there is some belief that the operator a may 'destruct' quanta (photons) in an arbitrary state \ifi). 
This belief is reflected even in the name 'photon-subtracted state', sometimes used for the state a\ip) |(| 0, ||. 
However, this concept, guided maybe more by intuition than by a sound proof, seems to be misleading when 
one deals not with single Fock states, but with their superpositions or quantum mixtures. 

The better known counter-example is the harmonic oscillator coherent state |a): since it is an eigenstate 
of a, the state a|a)/|a| has the same mean number of photons as |a) (see a similar remark in |Q). But as 
we will see below, this example is not the only one. In contradistinction, when operator a) is applied on 
an arbitrary state \ip) it produces a new one, whose mean number of quanta is always greater than that of 
therefore justifying its designation as 'creation operator'. Our aim is to provide a deeper investigation 
of the problem of 'photon subtraction' or 'photon addition' for arbitrary quantum states, drawing readers' 
attention to apparent 'paradoxes' and giving possible explanations. 

The plan of the paper is as follows. In Section 2 we discuss the relations between the 'photon excess' in 
the 'photon-subtracted' states and Mandel's parameter, illustrating them in the example of superposition of 
two Fock states. The concept of 'hyper-Poissonian' states is introduced in Section 3. The examples of such 
states include certain superpositions of the coherent and vacuum states and the negative binomial states. 
New families of 'logarithmic' states corresponding to the given values of the mean number of quanta and 
Mandel's parameter are considered in Section 4. A distinguished role of the exponential phase operators and 
their eigenstates - the phase coherent states - is discussed in Section 5. Section 6 concludes the paper. 
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2 'Photon-subtracted' states, 'photon excess' and Mandel's pa- 
rameter 



For any normalized state 

oo oo 

M=X>|n>, ]T|c n | 2 = l, (3) 

n=0 n=0 

the normalized 'photon-subtracted state' is given by 

_^ _^ oo ^ oo 

IV* — > = -j=a\i/>) = — == V Cn-s/ri \n - 1) = — = V c ri+ iV«+ (4) 
yn yn \ yn *—t 

n—l n—0 

where 

oo 

n = (^>|n|<0) = ^n|c„| 2 (5) 

n=l 

is the mean value of the quantum number operator h = a' a in . (In this article we consider mainly pure 
states, with the only aim to simplify notation.) Designating the mean number of quanta in the state \tp-) as 

oo 

N- = {rl>-\h\4>-) = - y>(n-l)|c„| 2 , (6) 

n=l 

it is easy to see that the difference between N- and n (the 'photon excess') is nothing but Mandel's q- 
parameter jjj: 

iV_ — n = n 2 /h — 1 — n = q . (7) 

Note that the transformation \tp) — > |^>_) is not a pure mathematical exercise. It is a real consequence 
of conditional measurements on beam splitters |(| [ic}| . Moreover, according to the Srinivas-Davies theory of 
photodetection |Tl| , the quantum states of field modes undergo the same transformation upon detection of 
one photon. In the last case, equation ([?]) and the related interpretation of the physical meaning of Mandel's 
parameter were discussed in ^| . 

As follows from (^), the 'annihilation' operator a diminishes effectively the mean number of quanta (q < 0) 
only for states having sub-Poissonian statistics (which have been always considered as 'nonclassical'). There 
are many families of 'sub-Poissonian' states, besides the Fock ones. In particular, all 'binomial' states jl4[ |l5| 
possess sub-Poissonian photon statistics. Other examples include the families of so-called 'Barut-Girardcllo 
states' or odd coherent states fl8|| . 

For any super-Poissonian state, the action of the 'annihilation' operator results in increase of the mean 
number of quanta. The explanation of such a 'paradoxical' behaviour stems from the first factor, \/n, in the 
right-hand side of equation ([!]): besides shifting the Fock state |n) by one quantum to the left, the operator 
& also increases the amplitude of the new state. Therefore, under certain conditions the relative weight of 
Fock states with larger numbers of quanta in state (^) can be greater than in the initial state ([|), which 
results finally in the increase of the mean number of quanta (cf . [ fL2| ) . 

2.1 Superposition of two Fock states 

Suppose the superposition of two different Fock states 



\iP)f =Vr\n) +Vl ~r\m), < r < 1. (8) 

Then 

, . „. rn(n — 1) + (1 — r)m(m — 1) 

n = rn + (1 - r)m, iV_ = — — y - -, 

rn + (1 — r)m 

and the condition iV_ > n is equivalent to the inequality r(l — r)(n — to) 2 > r(n — m) + m, which can be 
satisfied for any r (different from zero and I) if the difference n — m is sufficiently large. 

Since there are no limitations on the positive values of Mandel's parameter, applying the 'annihilation' 
operator to highly super-Poissonian states one can increase the mean number of quanta by any desired 
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quantity. For example, in the case of state (||) one has q f=a (1 — r)(n — m) if r(n — m) 3> m, and one can 
obtain g 3> 1 if also (1 — r)(n — m) 3> 1. On the other hand, by calculating the g-factor for the state \ip-) f, 
one can see that it tends to the limit value 

1 - r 

q°° = m- 1, 

r 

when n — > oo, for fixed m and r. So, if (1 — r)m <C 1, but (1 — r)n 3> 1, then by applying the operator a 
on one can significantly increase the mean number of quanta and simultaneously transform the highly 
super-Poissonian state into sub-Poissonian state \tp-}, which turns out to be much more closer to the 
state \n— 1) than the initial state was, with respect to \n). This happens due to a significant reduction of the 
relative weight of the low-energy state. Indeed, if initially this weight was equal to (1 — r) /r, in the new state 
it becomes approximately n/m 3> 1 times smaller. In other words, the operator a effectively 'annihilates' 
the low-energy components of quantum superpositions, increasing the relative weight of the higher-energy 
components. 



3 'Hyper-Poissonian' states 

Acting on the state \ip) by the operator a) we obtain the 'photon-added state' H ||, 

I 1 °° 

\^ + ) = n ==a^) = n ==Y,Cnsfl + ^\n + l). (9) 
v 1 + n v 1 + n Z~o 

The mean value of the quantum number operator in the state \ip+) equals 

N + = (V + |n|V + ) = f> + I) 2 M 2 = n + 1 + { -^. (10) 

1+71 1+71 
n— 

The difference N + — n is always not less than 1 (being equal to 1 only for the Fock states). Therefore the 
name 'creation operator' is justified for the operator a) . Nevertheless, in certain cases one can add much 
more quanta by acting on state \ip) not with the 'creation' operator, but with the 'annihilation' partner! 
Using the relation (An) 2 = 7l(l + q) together with equations (Q) and (|Io| ) one can see that the condition 
N- > iV+ is equivalent to 

q>l + 2n. (11) 



The states possessing the property q > n were named in 1 20 , |2lJ 'super-chaotic' or 'super-random' states 
(for such states, the second-order correlation parameter fir^(0) is greater than its value for the thermal 
states 5( 2 e rm (0) = 2). Therefore the states possessing the property ( |TT| ) can be named 'hyper-chaotic' or 
'hyper-Poissonian'. Note that Mandel's parameter of the squeezed vacuum states is given exactly by the 
right-hand side of equation ( |i"l| ) (it can be calculated, e.g., using the results of Consequently, the 

'hyper-Poissonian' states can be characterized as those for which the fluctuations of the photon number are 
stronger than in the squeezed vacuum states. 

3.1 Superposition of coherent and vacuum states 

Let us consider a superposition of the coherent state \a) and the vacuum state |0) 

m 0a = Vv\<x)+Z\0), 77+|£| 2 + 2^,e- Q2 / 2 Ree = l. (12) 

Since the overall phases are not essential, we assume that rj and a are real positive numbers (whereas the 
coefficient £ may be complex; however, it does not enter the formulae which we are interested in). In this case 
the operator a 'annihilates' the vacuum component of the initial superposition (p^), transforming \ip}o a into 
the coherent state \ip-)o a = \a). This transformation is accompanied by an increase of the mean number of 
quanta under certain conditions: 

n = 77a 2 , 7V_=a. N + = - — '— . (13) 

1 + na z 
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The ratio N + /N- equals 

N + _ j 1 - a- 2 - 3n 
N- ~ 1 + na 2 

thus for 3i] + a~ 2 < 1 we have N + < N_. Moreover, if a > 1, but r/a 2 <C 1, then n <C 1 < N + <C iV_. The 
probabilities of finding n quanta in the superposition (|l^) are 

1 - e~ a ) , p„ = r^e"" — . (14) 

A special case of the state ( [l2| ) without vacuum component (po = 0) was considered in |2j| as an example 
of 'the most nonclassical' state (see also a comment in 



3.2 Generating functions and negative binomial states 

It is well known that the statistical properties of the probability distribution {p n = \c n \ 2 } can be calculated 
with the aid of the generating function 



oo 

G(z) = y>„z", Pn =- j G^(0), G(l) = l, 



(15) 



since its derivatives at z = 1 generate the factorial moments, 



d r G{z)/dz r \ z=1 = n(n - l)(n - 2)...(n - r + 1) = n 



(r) 



(16) 



If the function G(z) for the state 1-0) is known, then the functions G±(z) for the states \ip±) are given by 
the formulae 

G-(z) = G + (z) = J— g ±[zG(z)]. 

n az 1 + n dz 

Applying these relations to the super-Poissonian negative binomial state |26| 



n=0 



1/2 



l«>, 



M > 0, <£ < 1, 



(17) 



(18) 



for which 



we find 



£/(!-£)) n = pq, 



G{z) 



G_(*) = 



1-g 



G+(z) = z 



1-C 



i + gjOf - 1) 

1 + ' 



Consequently, applying the operator a to the state ([^J) several times, one can add each time equal portions 
of quanta q (N- = fiq + q = h + q) without changing Mandel's parameter (g_ = q). On the other hand, 
after applying the operator a* on \tp > the mean number of quanta is increased, 



n = 1 + 



and one can easily verify that for fi < 1/2 < £(1 — /i) the action of the 'annihilation' operator adds more 
quanta than the action of the 'creation' operator. 



4 Quantum states with given mean photon number and Man- 
del's parameter 

Suppose that we want to find a probability distribution {p n } with two given parameters: the mean number 
n and the 'mean photon multiplication factor' 7 = N-/n. Since iV_ = alG- /dz\ z=1 , we can rewrite the 
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second condition, using the first of the relations in equation (|T7|) , as G"(l) = jhG'(l). Assuming that this 
relation must hold not only for z = 1, but also in some interval of values of the auxiliary variable z, including 
the point z = 1, we obtain a simple equation, G"(z) = 7nG'(z), whose solution, satisfying the conditions 
G(l) = 1 and G'(l) = n, reads 

G(z) = l- -+ -e^-V. (19) 

7 7 

The corresponding 'photon number' distribution is 

po = 1 - - + -e" 7f \ p„ = ~e-T» n > 1. (20) 

77 7 n! 

A pure quantum state possessing the photon number probability distribution {p n } can be written as 

oo 

M=y>*»v^i»>, (2i) 



n=0 

where the set of phases {</> n } may be quite arbitrary In particular, one can adjust the phases in such a 
way that function together with the weight set (20) become a special case of the superposition (|l2|), 



with |a| 2 = 7n and r\ = 7 _1 . The probability distribution ( p0| ) exists for any value of n if 7 > 1. For 
7 < 1 (sub-Poissonian case) it has sense only for not very large values of the initial mean photon number: 
jn < I ln(l — 7)|. For 7 — > 1 the Poisson distribution is recovered, while for 7 — ► 00, po = 1 and all other 
probabilities go to zero. 

Of course, there are infinitely many different distributions resulting in the same values of two parameters, 
n and 7 or n and q — (7 — l)n. For example, one can start from the equation G"{z) — jn 2 G(z), which 
results in the function 

G(z) = cosh [n^/7 (z - 1)] + 7~ 1/2 sinh [n^j (z - 1)] . (22) 

However, this function can be interpreted as generating function for nonnegative probabilities {p n } only 
under certain restrictions on the values of h and 7, because the coefficient pi in the Taylor expansion ( |l5|) of 
function (|22| ) becomes negative for n^j 3> 1 and 7 > 1, whereas for po the same happens if n^/7 3> 1 and 
7 < 1. This example shows the main difficulty of finding new distributions through differential equations 
for their generating functions: although one can write infinitely many differential equations for G{z) which 
result in the necessary conditions at z = 1, only few of them can be solved analytically, and even much 
fewer functions thus obtained have all derivatives at z = nonnegative. Nonetheless, several interesting 
distributions can be found in this way. 

4.1 Generalized logarithmic states 

For example, looking for probabilities resulting in a given value of the 'photon excess' (or Mandel's parameter) 
q 7 we can rewrite equation (Q) (or n 2 — n 2 — (q + l)n = 0) as (using relations jl^)) 

G"(l)- [G'(l)f -qG'(l) = 0. (23) 

One of many possible extensions of this equality to a finite interval of values of variable z is the equation 

G"{z) - [G'{z)] 2 - qG'(z) = 0. (24) 

Although being nonlinear this equation is easily solved by several methods, for example, multiplying it on 
the left by the integrating factor e~ G one gets a linear equation dF/dz = qF for F = e~ G G' . Taking into 
account the conditions G(l) = 1 and G'(l) = n, we obtain 



G(z) = 1 - In 

The probability distribution for q > is 

qe 



l + 2(l_ e 9(*-i) 



(25) 



Po = In 



q + n(l- e-i) 
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Although all coefficients p n with n > 1 are positive for any positive q and n, the vacuum probability po is 
nonnegative provided 

^<^. (27) 

1 — e 9 

For g >• 1 and q >• n the probabilities (p6|) become close to (pfj): 



1,2,3,. 



Using equations (|T7J) and ( f[7| ) we obtain Mandel's parameter q_ for the state |?/>_ 
factorial moments of the initial arbitrary state 



,.(3) 



7.( 2 )1 2 



q- = 



(28) 

in terms of the 
(29) 



The peculiar feature of the state generated by function ( pq ) is that for this state g_ = n, i.e., while for state 
n is the mean photon number and q is Mandel's parameter, for state \ip-} the mean photon number is 
n + q and Mandel's parameter is q_ = n. And the 'excess of photons' of state a\ip-} is n. 
Going to the limit q — > in (E5[) we obtain the generating function 



G{z) = 1 -ln[l -n(z - 1)] 



which yields the following probabilities and factorial moments (16) 

1 - In (1 + ft) , 



Pa 



1 

Pn = - 

n 



1 



n > 1, 



(r - 1)! n r 



1,2, 



(30) 



(31) 



(32) 



The probability distribution with the same factorial moments ([32|) was found in p7| practically by the same 
method. The difference is that the authors of |^7j calculated the generating function not for the probabilities, 
but for the factorial moments, having obtained, in particular, the function 



Qi(as) = Yj(-x) v {a ]v a v )/v\ = (A- l)" 1 In [(1 - A)nx + 1] + 1, 



(33) 



where the positive constant A was considered as a measure of the relative deviation from the Poissonian case 
(the article |27| appeared before H): 



[{a)a){a)a) - 1)) - {a)a) 2 } / \o) a) 2 = -A. 



(34) 



In the limit case A = formula (|33j) yields the factorial moments in the form ( p2|) (but for A ^ and 
q = —An 7^ the functions (^BJ) and (^) give different values of factorial moments). 

The statistics of the probability distribution (|l]) (not found explicitly in ^7j) is quite different from the 
Poissonian statistics, although this distribution gives Mandel's factor q = and the same first two factorial 
moments as in the coherent state with \a\ 2 = n. The corresponding quantum state was named in |^7j sub- 
coherent state (the same name was used in [[l7| for quite different states possessing sub- Poissonian statistics) . 
We can note in this connection that one of many possible pure states corresponding to the distribution ( |3l| ) 
is a special case of the family of so-called logarithmic states introduced in |2q] 



/log 



=1°) + E ns 

i V n 



= 1 + In (1 - 



1-|*| 5 



(35) 



G 



5 A distinguished role of the exponential phase operator and 
phase coherent states 

Acting on the 'logarithmic' state ( |35|) by the operator d we arrive at its 'subtracted' partner 

oo 

iWo^v^E^w- ( 36 ) 



n=0 



known under the names coherent phase state |30| , |3l[ , harmonious state P4I , or pseudothermal state f33| 
(see also @). 

The mean number of quanta in the 'sub-coherent' state ( |35| ) cannot exceed a small value e — 1 ~ 1.72. 
However, it is easy to show that there exist many sub-coherent states {q — 0) with arbitrary values of n. An 
example is the superposition of two Fock states (|J). Considering for simplicity the case of n ^> m, one can 
verify that q = either for r ~ m/n 2 (when n ps to) or for 1 — r 1/n (when n ~ rt). 

The coherent phase state ([56]) is the eigenstate of the exponential phase operator 

00 

E- = \ n - W = (a^)" 172 a = (n + 1)~ 1/2 a, (37) 

n=l 



introduced in |3^] and discussed in |3C], |3l], 36 1. Its Hermitian conjugated partner is 



B+ = £|n)<n-l| = 6 t (aa t ) 1/2 = d f (n + 1)" 1/2 , (38) 

n=l 

[B_,^ + ] = |o)<o|, Jsr+JS- = 1 - |o><o|. 

Acting on the Fock states, the operators E± shift the number of quanta by ±1 without changing the amplitude 
of the state vector: 

E_\n) = (l-6 n0 )\n-l), E + \n) = \n + l). (39) 

Therefore, applying the operators E± to an arbitrary state \ip) (|3|) we obtain the following normalized states 
$±)~E±\1>): 

-OO OO 

hM= 7 === Y,^n+M), $+)=Y i c n \n + l). (40) 

V 1 _ l C 0| n=0 n=0 

Their quantum number generating functions are related to the initial generating function G(z) ([l5]) by simple 
relations 

G-(z)= G{z) - G{ °\ G + (z) = zG(z). (41) 

Z(L-PO) 



Consequently, the mean numbers of quanta in the new states ( flpp are connected with the mean number n 
in the initial state as follows: 

N_ = (^_\h\if_) = - 1, N + = $ + \h$+)=n+l. (42) 
1 -po 

We see that the operator E+ adds exactly one quantum to any quantum state, whereas the operator _E_ 
removes exactly one quantum from any state which has no contribution of the vacuum state (i.e., if po = 0). 
Therefore, perhaps, not operators a and a', but operators £L and E + should be named 'annihilation' and 
'creation' operators in the literal meaning of these words! (See in this connection study |p7| , where it was 
shown how the states E^lip) could arise as a result of the interaction of the Jaynes-Cummings type between 
atoms and cavity fields or between motional and internal degrees of freedom of trapped ions. Shifted thermal 
states, which can be written as pf^ 1 ^^ — E™ pthE™ , have been considered in (3^], whereas methods of 
generating such states in a micromaser were discussed in (39).) 
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For the coherent states po = exp(— n), and for nClwe have iV_ = n/2. For the coherent phase states 
( |36| ) we have N- = h for any n = \z\ 2 / (l — \z\ 2 ). For many other states with po > 0, one can always make 
the 'modified photon excess' 

gr s #__ ft= (43) 
1 ~Po 

positive by increasing the value of n, which is, in a generic case, a free parameter independent of po- This 
can be easily seen, e.g., in the example (|2^), when the quantity 

(l-p )q = n- 2+1 (l- e -^) 
7 

is obviously positive for sufficiently large n and 7 > 1. The same is true for the negative binomial states 
( |l8| ) with jU < 1 and £ sufficiently close to 1. Moreover, if h > 2 (1 — po) /Po-, then N- > N + . 

Even using operators of the form A — f(h)a we can find the states \ip) for which state A\ip) has more 
quanta (in average) than the state |V>), for arbitrary function f(n). The simplest example is the 'two-Fock' 
state (Q) with m — 0. In this case n = rn. Any operator f{n)a transforms this superposition into the single 
Fock state \n — 1), and obviously (n — 1) > rn if (1 — r) > l/n. 



6 Conclusion 

Concluding, we have shown that the standard 'annihilation' boson operator diminishes the mean number of 
quanta only for 'nonclassical' states with negative values of Mandel's parameter, whereas there are many 
states (all of them superpositions or mixtures of the Fock states), with the peculiarity that when acted 
by a produce new states whose mean number of quanta is increased. Moreover, we have shown that in 
certain cases the 'annihilation' operator a can increase the mean number of quanta much more effectively 
than its 'creation' partner . We emphasized an interpretation (discovered in [fl2| , [l3f but still not well 
known) of Mandel's g-factor as the 'photon excess' in the mean number of quanta arising as a result of 
action of operator a on the given quantum state. Besides, we have shown that maybe the pair of Hermitian 
conjugated exponential phase operators could be the best candidates for being called 'annihilation-creation' 
operators, although the 'lowering' operator from this pair also cannot diminish the mean number of quanta 
by exactly one unit, for quantum states with nonzero vacuum component. This gives a hint that, perhaps, in 
the quantum photocount formalism of continuous measurements, as for instance, that proposed by Srinivas 
and Davies jll] one should use instead of a and the operators E- and E+, which would not violate the 
property V of the cited work (namely, the condition of boundness for the counting rate). However, this issue 
is outside the objective of this letter, and we shall report on it elsewhere. 

SSM acknowledges partial financial support by CNPq (Brasilia, DF) and FAPESP (Sao Paulo; Contract 
No. 2000/15084-5). WD thanks CNPq for the full financial support. 
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